Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



The Cluster Substructure - Alignment Connection 



Manolis Plionis^ & Spyros Basilakos^ 

^Institute of Astronomy & Astrophysics, National Observatory of Athens, I.Metaxa & B.Pavlou, Palaia Penteli, Athens 152 36, Greece 
^Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Gonsort Road, London SW7 2BW, UK 



1 February 2008 



o 
o 

> 

o 

:z; 

o 

(N 
> 
m 



Oh; 

6 



1 INTRODUCTION 



An interesting observable, that was thought initially to pro- 
vide strong constraints on theories of galaxy formation, is 
the tendency of clusters to be aligned with their nearest 
neighbour as well as with other clusters that reside in the 
same supercluster (cf. Binggeli 1981; West 1989; Plionis 
1994; Chambers, Melott & MiUer 2001). Analytical and nu- 
merical work have shown that such alignments, expected 
naturally to occur in "top-down" scenarios (cf. Zeldovich 
1970), are also found in hierarchical clustering models of 
structure formation like the CDM (Bond 1986; West et al. 
1991; Splinter et al. 1997; Onuora & Thomas 2000). This 
fact could be explained as the result of an interesting prop- 
erty of Gaussian random fields that occurs for a wide range 
of initial conditions and which is the "cross-talk" between 
density fluctuations on different scales. Furthermore, there is 
strong evidence that the brightest galaxy (BCGs) in clusters 
is aligned with the orientation of its parent cluster and even 
with the orientation of the large-scale filamentary structure 
within which they are embedded (cf. Struble 1990; West 
1994, Fuller, West & Bridges 1999). 

Within the framework of hierarchical clustering, the 
anisotropic merger scenario of West (1994), in which clus- 
ters form by accreting material along the filamentary struc- 
ture within which they are embedded, provides an inter- 
esting explanation of such alignments as well as of the ob- 
served strong alignment of BCGs with their parent cluster 
orientation. In this framework, one should expect that dy- 
namical young clusters, at there early stages of formation 
in which they are not smooth and spherically symmetric, 
should show substructures that are aligned with the local 
large-scale structures, an effect observed also in numerical 
simulations for a variety of power-spectra (van Haarlem & 
van de Weygaert 1993; Tormen 1997). Indeed, supporting 
this view West, Jones & Forman (1995) found, using 43 
Einstein clusters that have a neighbour within ^ 10 h'^ 
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Mpc, that cluster substructures do show a tendency to be 
aligned with the orientation of the major axis of their parent 
cluster and with the nearest-neighbouring cluster (see also 
Novikov et al 1999). In this work we investigate the relation 
between the strength of cluster-cluster alignments and the 
large-scale environment in which the clusters are embedded 
using the APM cluster catalogue (Dalton et al 1997), which 
is the largest one available. 



2 METHODOLOGY 

The APM cluster catalogue is based on the APM galaxy 
survey which covers an area of 4300 square degrees in the 
southern sky containing about 2.5 million galaxies brighter 
than a magnitude limit of hj = 20.5 (for details see Mad- 
dox et al. 1990). Dalton et al (1997) applied an object cluster 
finding algorithm to the APM galaxy data using a search ra- 
dius of 0.75 Mpc in order to minimize projection effects, 
and so produced a list of 957 clusters with Zest^ 0.13. Out 
of these 309 (~ 32%) are ACQ clusters, while 374 (~ 39%) 
have measured redshifts (179 of these are ACQ clusters). 
The APM clusters that are not in the AGO list are rela- 
tively poorer systems than the Abell clusters, as we have 
verified comparing their APM richness's (see Dalton et al 
1997 for definition of richness). 

For the present analysis we use 903 of the above APM 
clusters, since 54 clusters are found in the vicinity of plate- 
holes or crowded regions, a fact which affects severely their 
shape parameters. The cluster distance is estimated from 
their redshift using Ho = lOO/i km s^^ Mpc"^ and qo — 0.5. 



2.1 Cluster Shape Parameters and Alignment 
Measure 

A detailed analysis of the cluster shape determination pro- 
cedure and of the intrinsic APM cluster shapes can be 
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Figure 1. Distribution of APM cluster position angles. 



found in Basilakos, Plionis & Maddox (2000). Here we only 
sketch the basic procedure which is based on the familiar 
moments of inertia method with Jn — ^ Wi(r1 — x1), 
h2 = Wii^rf - yf), Ii2 = hi = - J2 ^i^ilH^ where 
Xi and yi are the Cartesian coordinates of the galaxies that 
their projected separation is such that they are judged as 
belonging to the cluster (details in Basilakos et al 2000) and 
Wi is their weight. We, then diagonalize the inertia tensor 
solving the basic equation: 



det{hj ~ y Mi) = 0, 



(1) 



where M2 is the 2x2 unit matrix. The cluster ellipticity is 
given by e = 1 — where Xi are the positive eigenvalues 
with (Al > A2). This method can be applied to the data us- 
ing either the discrete or smoothed distribution of galaxies. 
The determination of the cluster orientation is consistent 
among the two methods but this is not true also for the 
cluster ellipticity (for details see Basilakos et al. 2000) 

In Figure 1 we present the derived cluster position an- 
gle, 9i, distribution. We see no significant systematic ori- 
entation effects, a fact that we quantify estimating the 
Fourier transform of the galaxy position angles; Cn = 
^ cos 2n9i and Sn = \/^ X/ 2n9i . If the galaxy 
position angles are uniformly distributed between 0° and 
180°, then both C„ and Sn have zero mean and unit stan- 
dard deviation. Therefore large values (> 2.5 — 3) indicate 
significant deviation from isotropy. We find that Ci,2 and 
Si, 2 have values ^ 1.2, which indicates that there is no sys- 
tematic orientation bias. 

In order to investigate the alignment between cluster 
orientations, we define the relative position angle between 
the major axis orientation of a cluster and the direction 
to a neighbouring one by = l^i — (where is the 
position angle of the cluster pair separation vector). In an 
isotropic distribution we will have {(j>i,j) — 45°. A significant 
deviation from this would be an indication of an ainsotropic 
distribution which can be quantified by (Struble & Peebles 
1995): 



N 



45 



(2) 



0, while the stan- 



In an isotropic distribution we have (5) 
dard deviation is given by a = 90/vl2A'^. A significantly 
negative value of S would indicate alignment and a positive 



misalignment. To avoid problems related to ill defined posi- 
tion angles we will use in our analysis only clusters that have 
ellipticities > 0.05 (ie., iV = 888)Q It should be noted that 
systematic biases and projection effects will tend to mask 
any true alignment signal. For example, the projection of 
foreground galaxies along the line of sight of a cluster as well 
as the projection on the plane of the sky of member galaxies, 
always work in the direction of smearing alignments. 



2.2 Substructure Measure &: Significance 

Major obstacles in attempting to determine the dynami- 
cal state of a cluster is (1) the ambiguity in identifying 
cluster substructure in 2D or even 3D cluster data and 
(2) the uncertainty of post-merging relaxation timescales. 
Evrard et al. (1993) and Mohr et al. (1995) have sug- 
gested as an efficient indicator of cluster substructure the 
shift of the center-of-mass position as a function of den- 
sity threshold above which it is estimated. The centroid- 
shift (sc) is defined as the distance between the cluster 
center-of-mass, {xo,yo), which may change at different den- 
sity thresholds and the highest cluster density-peak, (a::p, j/p), 
ie., sc = ^y {xo — Xp)^ + (xo — Xp)^. 

Kolokotronis et al. (2001), using in a complementary 
fashion optical and X-ray data (see also Rizza et al. 1998), 
since in the X-ray band projection effects are minimal, cali- 
brated various substructure measures using APM data and 
pointed ROSAT observations of 22 Abell clusters and found 
that in most cases using X-ray or optical data one can iden- 
tify substructure unambiguously. Only in ~ 20% of the clus- 
ters that they studied did they find projection effects in the 
optical that altered the X-ray definition of substructure. An 
important conclusion of Kolokotronis et al. (2001) was that 
a large and significant value of sc is a clear indication of 
substructure in APM optical cluster data. 

The significance of such centroid variations to the pres- 
ence of background contamination and random density fluc- 
tuations are quantified using Monte Carlo cluster simula- 
tions in which, by construction, there is no substructure. 
For each APM cluster, 1000 simulated clusters are pro- 
duced having the observed ellipticity, the observed number 
of galaxies, following a King's profile, as well as a random 
distribution of expected background galaxies, determined by 
the distance of the cluster and the APM selection function 
(note that the number of "galaxies" that we allow to follow 
the King's profile is the observed number minus the expected 
random background). The King-like profile is: 



E(r-) oc 



1 



{0 



(3) 



where rc is the core radius. We use the weighted, by the sam- 
ple size, mean of most recent rc and a determinations (cf. 
Girardi et al. 1998), i.e., rc ~ 0.085 fe"^ Mpc and a ~ 0.7. 
We do test the robustness of our results for a plausible range 
of these parameters (details can be found in Kolokotronis et 
al. 2001). 

Naturally, we expect the simulated clusters to generate 
small sc's and in any case insignificant shifts. Therefore, 



* Our results remain unaltered even for higher values of the el- 
lipticity cutoff. 
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from each set of Monte-Carlo cluster simulations we derive 
{sc)sim as a function of the same density thresholds as in the 
real cluster case. Then, within a search radius of 0.75 
Mpc from the simulated highest cluster peak, we calculate 
the quantity: 

^ _ {sc)o — {sc)sim j.^^ 

which is a measure of the significance of real centroid shifts 
as compared to the simulated, substructure- free clusters. 
Note that (sc)o is the average, over three density thresh- 
olds, centroid shift for the real APM cluster. 

A further possible substructure identification procedure 
is based on a friend-of-friends algorithm, applied on 3 over- 
density thresholds of each cluster (for details see Kolokotro- 
nis et al. 2001). Three categories are identified, based on the 
subgroup multiplicity and size: (a) No substructure (uni- 
modal), (b) Weak substructure (multipole groups but with 
total group mass < 25% of main), (c) Strong substructure 
(multipole groups but with mass > 25% of main). 



3 RESULTS & DISCUSSION 

3.1 Cluster Substructure and Alignments 

Applying the centroid shift substructure identification proce- 
dure to the 903 APM clusters we find that about 30% of clus- 
ters have significant (> 3cr) substructure. Note that defin- 
ing as having significant substructure those clusters with 
a > 2.5 or 2 increases the fraction to ~ 40% and 50% re- 
spectively. Furthermore, changing the structural parameters 
of the Monte-Carlo clusters changes the actual cr-values, al- 
though their relative significance rank-order remains unal- 
tered. Alternatively if we apply the subgroup categorization 
procedure we find that ~53% of the APM clusters show 
strong indications of substructure. 

We have tested whether there is any systematic red- 
shift dependent effect of the cluster substructure categoriza- 
tion and found none (using either measured or estimated 
redshifts). We conclude that irrespectively of the method 
~ 30% — 50% of the APM clusters show indications of sig- 
nificant substructure (in accordance with Kolokotronis et al 
2001). 

We now test whether the well known nearest-neighbour 
alignment effect, present in the Abell clusters (cf. Bingelli 
1982; Plionis 1994), is evident also in the poorer APM clus- 
ters. In table 1 we present, for the whole APM cluster sam- 
ple, our alignment results as a function of maximum inter- 
cluster separation, Dmax, for both the nearest-neighbours 
(Snn) and for all neighbouring pairs (San), within Dmax- 

It is evident that there is significant indication of cluster 
alignments, with the alignment signal dropping in amplitude 
and significance, as a function of increasing Dmax- In order 
to test whether this result is dominated by the ACO cluster 
pairs, and thus whether it is a manifestation of the already 
known Abell cluster alignment effect, we have excluded such 
pairs (113/888 for the nn-case) to find consistent but more 
significant alignment results. For example for the Dmax = 15 

Mpc case, we obtain (Snn) — — 4.3±1.3 with P{> x^) = 
0.0004 and (San) ^ -2.7 ± 0.9 with P(> x^) = 0.001. 

Note that a large number of APM clusters have esti- 
mated redshifts and thus the cluster distance uncertainties 



Dmax 


(Snn) 




^ pairs 


(San) 




^ pairs 


5 


-7.7±2.7 


0.005 


90 


-7.5 ±2.5 


0.002 


110 


10 


-4.6±1.6 


0.005 


270 


-4.1±1.3 


0.000 


374 


15 


-3.8±1.2 


0.001 


444 


-2.6±0.9 


0.002 


848 


20 


-2.9±1.0 


0.004 


616 


-1.1±0.7 


0.023 


1602 


30 


-1.7±0.9 


0.069 


805 


-0.9±0.4 


0.112 


4012 



Table 1. Nearest neighbour (rm) and all-neighbour (an) align- 
ment signal as a function of maximum pair separation, Dmax- 
The probabilities are derived by comparing the binned 0- 
distribution (using 3 bins of 30° width) with the Poisson expec- 
tation values. 



N . N . 

pairs pairs 

1602 916 540 278 164 58 40 10 801 458 270 139 82 29 20 5 

' I ' ' ' I ' ' ' I ' ' 




4 6 2 4 

a a 

Figure 2. Alignment signal of all cluster pairs with separation 
Dec < 20 h^^ Mpc, as a function of substructure significance, 
a- The left panel presents alignments between cluster major axis 
orientation and the direction to a neighbour while the right panel 
represents the alignments between the two cluster major axes 
orientation. Furthermore, the filled symbols represent the signal 
based only on the centroid- shift substructure categorization while 
the open symbols represent the signal from clusters that are also 
categorized as having strong substructure by the subgroup cate- 
gorization procedure. 



will tend to hide true alignments. Therefore, the measured 
alignment signal should be considered rather as a lower limit 
to the true one. 



3.2 The Substructure- Alignment Connection 

We have correlated the alignment signal with the substruc- 
ture significance indication in order to see whether there is 
any relation between the large-scale environment, in which 
the cluster distribution is embedded, and the internal cluster 
dynamics. 

In figure 2 we present the alignment signal, (5) , between 
all cluster pairs with separations < 20 Mpc that have 
substructure significance above the indicated a value. Evi- 
dently, there is a strong correlation between the strength of 
the alignment signal and the substructure significance level 
(the two panels are based on slightly different definition of 
the alignment signal - see caption for details). 
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r (h-' Mpc) 

Figure 3. Two-point correlation function of all APM clusters 

(open symbols) and of cr > 4 clusters (filled symbols). The lines 
represent the best (r /ro)~^'^ fit with ro — 12 and ~ 20 Mpc 
respectively. Insert: The cluster correlation length as a function 
of substructure significance. 

In order to assess the statistical significance of obtaining 
such & 5 — a relation, we perform a large number (10000) of 
Monte-Carlo simulations in which we reshuffle the measured 
position angles, assigning them at random to the clusters. 
Then we derive the corresponding 5 — a relation and count in 
how many such simulations do we get values of 5 which axe 
as negative, or more, than the corresponding ones of figure 2. 
We find that the corresponding probability is < 10~^ , while 
for the less restrictive case, of having negative (5-values (of 
any amplitude) for all cr, the probability is again quite small 
(~ 7.5 X 10"^). 

Note that from the analysis of Kolokotronis et al (2001) 
it is expected that our procedure will misidentify the dy- 
namical state of ~ 20% of the APM clusters. However, such 
misidentification will act as a noise factor and will tend 
to smear any true alignment-substructure correlation, since 
there is no physical reason why random projection effects, 
within 0.75 h^^ Mpc of the cluster core, should be corre- 
lated with the direction of neighbours within distances up 
to a few tens of Mpc's (such a correlation could be expected 
at some level only for nearest-neighbours in angular space 
but wc have verified that choosing such pairs wc obtain an 
insignificant alignment signal). Therefore our result, based 
on the largest cluster sample available, supports the hierar- 
chical clustering scenario and in particular the formation of 
cluster by anisotropic merging along the filamentary struc- 
ture within which they are embedded (cf. West 1994; West, 
Jones & Forman 1995). 

3.3 Local density - Substructure Correlation 

If the above view is correct then one would expect that clus- 
ters with significant substructure should be residing prefer- 



entially in high-density environments (superclusters), and 
this would then have an imprint in their spatial two-point 
correlation function. 

In figure 3 we present the spatial 2-point correlation 
function of all APM clusters (open symbols) and of those 
with substructure significance cr > 4 (filled symbols). It is 
clear that the latter are significantly more clustered. This 
can be seen also in the insert of figure 3 were we plot the 
correlation length, ro, as a function of cr, which is clearly 
an increasing function of cluster substructure significance 
level. Note that in order to take into account the possible 
systematic distance dependent effects in the diflFerent cluster 
subsamplcs wc generate random catalogues, used to normal- 
ize the number of cluster-cluster pairs in the ^(r) estimate, 
using the individual distance distribution of each subsample 
and not the overall APM cluster selection function. Had we 
used the latter, the increase of ro with a would have been 
more severe. 

Furthermore, we have tested whether this effect could 
be due to the well-known richness dependence of the cor- 
relation strength, and found a weaJj, if any, such richness 
trend. In order to further investigate this issue we have cor- 
related the APM riclmcss (sec Dalton ct al 1997) with a 
and found a very weak but significant correlation (the Pear- 
son's coefficient is only 0.26 but with a significance >99.9%). 
Therefore, a component contributing to the increase of the 
ro(cr) function could possibly be the cluster richness effect 
(through the above richness-a weak correlation). However, 
the increase of ro(cr) is quite dramatic and cannot be at- 
tributed only to this weak richness dependence. This is cor- 
roborated also from the fact that there is a larger fraction 
of clusters with significant substructure (high o values) re- 
siding in high-density regions (superclusters), as shown in 
figure 4 where we plot the fraction of clusters with cr > 3 
that belong to superclusters identified by the indicated per- 
colation radius. Evidently the fraction increases inversely 
with percolation radius. These results are similar to those of 
Schombert & West (1990) based on Abell clusters, in which 
they found that flattened clusters, flatness being an indicar 
tion of dynamical youth, are more frequent in high density 
environments (see however Herrera & Sanroma 1997 for a 
different view). 

The conclusion of this correlation function analysis is 
that indeed the clusters showing evidence of dynamical ac- 
tivity reside in high-density environments, as anticipated 
from the alignment analysis. It is interesting that such en- 
vironmental dependence has also been found in a similar 
study of the BCS and REFLEX clusters (Schiiecker et al. 
2001) and for the cooling flow clusters with high mass ac- 
cretion rates (Loken, Melott & Miller 1999). 



4 CONCLUSIONS 

We have presented evidence, based on the largest available 
cluster sample, the APM, that there is a strong link between 
the dynamical state of clusters and their large-scale environ- 
ment. Clusters showing evidence of dynamical activity are 
significantly more aligned with their nearest neighbours and 
they are also much more spatially clustered. This supports 
the hierarchical clustering models in which clusters develop 
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by accreting matter along the large-scale filamentary struc- 
tures within which they are embedded. 
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